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ABSTRACT 

A basic property of objects, like galaxies and halos that form in cosmological structure formation 
simulations, is their shape. Here, we critically investigate shape determination methods that are 
commonly used in the literature. It is found that using an enclosed integration volume and weight 
factors r~^ and r^j (elliptical radius) for the contribution of each particle or volume element in the 
shape tensor leads to biased axis ratios and smoothing of details when calculating the local shape 
as a function of distance from the center. To determine the local shape of matter distributions as 
a function of distance for well resolved objects (typically more than ©(lO"*) particles), we advocate 
a method that (1) uses an ellipsoidal shell (homoeoid) as an integration volume without any weight 
factors in the shape tensor and (2) removes subhalos. 
Subject headings: methods: data analysis — methods: numerical 
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1. INTRODUCTION 

Typically, the distribution of matter in objects that 
form in cosmological structure formation simulations is 
crudely described by spherically averaged density pro- 
files (e.g. Navarro et al. 1996; Moore et al. 1998). But 
real halos are not spherically symmetric and a natural 
extension is to describe the iso-density contours as sur- 
faces of ellipsoids. There is a wealth of literature with 
many different methods that are used to measure the lo- 
cal shape of a mass distribution (e.g. Gerhard 1983; Frenk 
et al. 1988; Katz 1991; Dubinski & Carlberg 1991; War- 
ren et al. 1992; Cole & Laccy 1996; Jing & Suto 2002; 
Springel et al. 2004; Kazantzidis et al. 2004; Bailin & 
Stcinmetz 2004, 2005; AUgood et al. 2006; Kuhlcn et al. 
2007; Bctt et al. 2007; Hayashi et al. 2007; Debattista 
et al. 2008; Warnick et al. 2008; Kazantzidis ct al. 2010; 
Kncbe ct al. 2010; Lau et al. 2011; Vera-Ciro et al. 2011). 
Their common goal is to recover the iso-density surfaces 
of the underlying matter distribution. Other character- 
istics, such as the potential, can also be used to describe 
the objects (e.g. Springel et al. 2004; Hayashi et al. 2007; 
Kazantzidis et al. 2010; Lau et al. 2011). 

Unfortunately, the literature lacks a systematic com- 
parison of the different methods - especially under con- 
trolled conditions where the exact shape is known. For 
some notable exceptions, see e.g. AUgood et al. (2006); 
Kazantzidis et al. (2004). But as far as we know, there 
is no publication that investigates the different methods 
under controlled conditions with known shape as it is 
done in this paper. Presumably, many of the quantita- 
tive discrepancies in the literature originate in the var- 
ious methods that are used for determining the shape. 
This work is intended to shed some light on the effects 
and systematics of the various methods that are based on 
an iterative procedure that uses a shape tensor with dif- 
ferent weighting schemes and integration volumes. The 
influence of the local mass density profile on the capabil- 
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ity of the shape finding method to recover the iso-density 
contours is also investigated. 

2. BACKGROUND 

From classical mechanics, the relation between the an- 
gular momentum vector L and the angular velocity vec- 
tor a; of a body is given by 

L=lLJ, (1) 

where I is the moment of inertia tensor defined by 



pir){r^l -rr^)dU , 



(2) 



where the integration is over the whole volume of the 
body and 1 is the identity tensor. Here, p{r) is the mass 
density at the location of the volume element dV pointed 
by the position vector r with respect to the center of the 
mass distribution. By defining the tensor 



M = I p{r)rr'^dV , 
Jv 



(3) 



which is the second moment of the mass distribution, it 
follows that 

/ = tr(M)i-M. (4) 

Hence, the tensor M is the fundamental quantity that 
describes how the matter is distributed. 
We now define the shape tensor as 

M Jyp{r)rr'^dV 



where 



Mtot IvPi^W 

Mt„t, = / p{r)dV 



(5) 



(6) 



is the total mass of the body. The shape tensor has units 
of length squared. For a discrete set of particles with 



Pi"^) = ^mfc(5(r- Tfc) 



(7) 
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Table 1 


Summary 


of methods 


Method 


w{r) 


V 


SI 


1 


eUipsoidal shell 


S2 


r 


ellipsoidal shell 


S3 


'^cU 


ellipsoidal shell 


El 


1 


enclosed ellipsoid 


E2 


r 


enclosed ellipsoid 


E3 


'■ell 


enclosed ellipsoid 



we obtain for the individual elements of the shape tensor 



S, 



Efc "^fc 



(8) 



where {rk)j denotes the j component of the position vec- 
tor of the fc-th particle and the summation is over all 
particles within the integration volume V. The tensors 
S and M describe how the mass is distributed, hence 
our choice for naming S the shape tensor. 

The tensors / and M have the same eigenvectors, lim 
is an eigenvalue of M, then tr(M') — m is an eigenvalue 
of I. The detailed meaning of the eigenvalues depends 
on the integration volume and the mass distribution (i.e. 
density profile). For example, for a thin ellipsoidal shell 
(a thin homoeoid) of uniform density, the eigenvalues 

of M are Meso.'^/S, Mes^V^ and MescVS (A^es is the 
mass in the ellipsoidal shell). Whereas for an ellipsoid 
of uniform density the eigenvalues are Meo^/S, Me^^^/S 
and Mec^/5 (Me is the mass in the elhpsoid). 

Unfortunately, the tensor M (Equation (3)) is often 
inaccurately denoted as the moment of inertia tensor in 
the astronomy and astrophysics literature. This proba- 
bly goes back to Binncy & Trcmaine (1987) (Page 494, 
Equation 8-11), where M was called the moment of iner- 
tia tensor. Fortunately, this was corrected in the second 
edition (Binney & Trcmaine 2008, page 796, Equation 
D-39). 

3. METHODS 

The shape tensor can be generalized by using an addi- 
tional weight function w{r) 



S ^ 



Jy p{r)w(r)rr^dV 
IvPir)dV 



(9) 



By setting w{r) — 1 and choosing p{r) to be the mass 
density we obtain our standard definition (Equation (5)). 
Other choices are also possible. For example, a weighting 
by number with p{r) = ^j, 5{r — r^) being the number 
density (which is, of course, equivalent to the mass den- 
sity weighting if all the particles have equal mass). Or 
p{r) — J2k '^('" ~ fk)/pk where pt is the local density of 
the particle like in Warnick et al. (2008). If one is in- 
terested in the shape of a matter distribution where the 
particles or volume elements can have a different mass 
(e.g. for gas and stars), it is essential to use p{r) as the 
mass density. Here, we only use p{r) as the mass den- 
sity. Throughout the paper, we use the elliptical radius 
Toll for distances from the center for ellipsoidal shapes. 
The elliptical radius r^w (see also Equation (10)) is the 
semi-major axis of the local homoeoid or ellipsoid. 



We concentrate on 6 different methods for determin- 
ing the shape of a matter distribution (see also Table 1). 
These methods differ by using a different integration vol- 
ume V and different weight functions w{r). For calculat- 
ing the local shape at a distance r^i, in the methods with 
a starting letter S, the integration is over an ellipsoidal 
shell (homoeoid) volume centered at reii (in logarithmic 
space). In the methods with first letter E, the integra- 
tion is over the whole enclosed ellipsoidal volume within 
Tell. For the different weight functions w{r), we use (1) 
w{r) — 1, (2) w{r) = r^^ and (3) w{r) = r, 
elliptical radius is given by 



2 „„A f'i\ .„t^\ - .-2^ The 



roll 




(10) 



where (aicii, j/oii, Zdi) are the coordinates of the volume el- 
ement or particle in the eigenvector coordinate system of 
the ellipsoid, i.e. Vgu corresponds to the semi-major axis 
a of the ellipsoid surface through that particle or volume 
element. Additionally, we also check for the importance 
of the removal of subhalos. Cases where we removed the 
subhalos are marked with a -, cases where they remained 
by a +. 

In order to calculate the local shape at a distance reii 
from the center, we use an iteration method (e.g. Katz 
1991; Dubinski & Carlberg 1991; Warren et al. 1992) and 
start with a spherically symmetric integration volume 
(shell or sphere). Then the shape tensor is calculated 
according to the different methods. By diagonalizing S 
we get the eigenvectors and eigenvalues at distance rdi. 
The eigenvectors give the directions of the semi-principal 
axes. The eigenvalues of S for the method SI are a^/3, 
6^/3 and c^/3 where a, b and c are the semi-principal 
axes with a > 5 > c - at least in the thin homoeoid 
approximation where the density is uniform. Hence, the 
square roots of the eigenvalues are proportional to the 
lengths of the semi-principal axes for method SI and we 
can readily calculate the axis ratios b/a and c/a. For 
method S3 we expect to get the same axis ratios as for 
method SI since dividing by the semi-major axis a — r^n 
squared, which is a constant for a thin ellipsoidal shell, 
just changes the geometrical meaning and normalization 
of the eigenvalues but not the axis ratios. 

For the other methods it is not clear what the detailed 
geometrical meaning of the eigenvalues is. For the meth- 
ods that use the enclosed ellipsoidal volume, this will also 
depend on the mass density profile. The r~^ weighting 
projects the volume elements dV onto the unit sphere. 
This projection complicates the physical interpretation 
of this method. It is generally assumed though that the 
eigenvalues of S in these cases are still proportional to the 
semi-major axes squared. Hence, we calculate the axis 
ratio for the other methods the same way as for methods 
SI and S3 - as it is generally done in the literature. 

We then keep the length of the semi-major axis fixed 
(but the orientation can change) and calculate S again by 
summing over all particles within the new deformed inte- 
gration volume (homoeoid or ellipsoid) with semi-major 
axis a = rdi and axis ratios b/a and c/a but with the 
new orientation. For the shape determination we allow 
volume elements or particles to be in several bins/shells. 
Of course this is naturally the case when using an en- 
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closed elhpsoidal volume. It is also necessary when using 
an ellipsoidal shell since neighboring shells can overlap 
due to slightly different orientation and axis ratios. This 
iteration is repeated until convergence is reached. As a 
convergence criterion we require that the fractional dif- 
ference between two iteration steps in both axis ratios is 
smaller than 10^'^. 

For methods using the shape tensor, it is important 
to use an iteration method that allows the algorithm to 
adapt the integration volume to the a priori unknown 
shape of the object. Often one also finds in the litera- 
ture that no iteration procedure is used and just a simple 
spherical shell or enclosed sphere is used as the integra- 
tion volume in order to calculate the shape (e.g. Ger- 
hard 1983; Frenk ct al. 1988; Cole & Lacey 1996; Bailin 
& Steinmetz 2004, 2005; Knebe et al. 2010). To us the 
physical meaning of the outcome of such a procedure is 
unclear and we do not further pursue it here. 

A further method for calculating the shape of contours 
is by selecting particles by their local density (e.g. Jing & 
Suto 2002; Warnick et al. 2008; Vera-Ciro et al. 2011) or 
potential (e.g. Springcl et al. 2004; Hayashi ct al. 2007; 
Kazantzidis ct al. 2010). There, no iteration procedure 
is needed. 

Often one also finds in the literature, that the moment 
of inertia tensor I (Equation (2)) in combination with an 
enclosed ellipsoidal integration volume is used for calcu- 
lating the axis ratios. This procedure assumes relations 
between the eigenvalues and semi-principle axes that are 
strictly valid only for a uniform ellipsoid or homoeoid 
(e.g. Bett et al. 2007; Warnick et al. 2008). For a thin 
homoeoid this is fine (under that assumption that the lo- 
cal density is constant in the shell) but for the enclosed 
ellipsoidal integration volume, the result is made equiv- 
alent to the method that just uses the shape tensor by 
construction. 

4. CONTROLLED CONDITIONS 

First, we examine the behavior of the different methods 
under controlled conditions where we know the correct 
shape. For this purpose, we set up various model halos 
that have different density, shape and orientation profiles 
with HALOGEN (Zcmp ct al. 2008). 

4.1. Models 

HALOGEN can generate random realizations of spheri- 
cal halos with a/37-profiles (Zhao 1996) 



p{r) 



Po 



(r/r,)'r[l-h(r/r,)"][(/^-^)/° 



(11) 



in equilibrium, where an importance sampling method 
(multimass technique) can be applied. Here, we just 
interpret the spherical radius r in the a/37-profiles as 
semi-major axis a = feii of a surface of an ellipsoid. For 
the generation of a uniform distribution of points on a 
surface of an arbitrary shaped ellipsoid, which is needed 
for setting up an ellipsoid with a given density profile, 
a method as outlined in Section 2.5.5 of Rubinstein & 
Kroese (2007) is used. Since we only care about the 
spatial distribution of the matter for our purpose, no ve- 
locities are assigned to the sampled particles. 
For the variation of the axis ratios and orientation with 
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Figure 1. Measured axis ratio b/a as a function of distance for 
halos with different inner slope 7 = (top panel), 1 (middle panel) 
and 2 (bottom panel). The halos were initialized with a constant 
axis ratio of b/a = 0.8 (thin dashed line). Except for the methods 
S2 and E2, all the other methods find the expected value. The 
fluctuations in the center seen for the inner slopes 7 = and 1 are 
mainly due to resolution and decrease when increasing the sampling 
resolution and using a finer binning. 



distance, we use a simple parametrization of the form 



x{r, 



clU 



loglo(^cll/?'0,2;) +2^0 



(12) 



Here, x can be b/a, c/a, 61, 62 and ^3, respectively. The 
angles 6*1, 6*2 and 6*3 are the Euler angles of an active 
z—x'—z" rotation. This allows us to twist the orientation 
of the principal axes as a function of distance. 

We use a generalized NFW (Navarro ct al. 1996) form 
for the density profile of the halos, i.e. we set a = 1, 
/3 = 3, and use 3 different values for the inner slope, i.e. 
7 = 0, 1 and 2. The ellipsoidal halos are sampled with 
10^ particles of the same mass within 10 r^ (no multi- 
mass technique applied). For some cases also different 
resolution halos with up to 10® particles within 10 r^ 
are used. To compare to current state-of-the-art cosmo- 
logical structure formation simulations: hydrodynamical 
simulations have reached 0(10^) particles per halo (e.g. 
Gucdcs et al. 2011; Zcmp et al. 2011) whereas halos in 
dissipationless N-body simulations are even resolved with 
O(IO^) particles (e.g. Springcl et al. 2008; Stadel et al. 
2009). Beyond 10 r^ an exponential cut-off of the mass 
density profile is applied in order to keep the total mass 
finite (for more details see Zcmp ct al. (2008)). With 
a resolution of 10^ particles, one can roughly sample an 
NFW profile down to 0.1 r^. The resolved scale depends 
on the inner slope 7. For 7 = 2, this scale is smaller and 
for 7 = it is larger (for more details see Zcmp et al. 
(2008)). Hence, for all profiles in the following plots only 
the range 0.1-10 r^ is shown. 

4.2. Constant axis ratios - aligned orientation 
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As a first deviation from perfect spherical symmetry, 
we set up halos with constant axis ratios, while the prin- 
cipal axes are kept aligned at all distances. Figure 1 
shows how the 6 different methods described in Section 
3 perform for our 3 halos with 7 = 0, 1 and 2. In these 
models we set b/a — 0.8 and c/a ~ 0.6. For clarity we 
only show the results for b/a. The findings are similar 
for c/a. 

The results for the methods SI, S3, El and E3 agree 
very well with the expected value. The small fluctuations 
seen in the center are due to resolution and depend on the 
mass profile as well. The fluctuations get smaller when 
sampling the same halo with more particles and using a 
finer binning. The default binning used in this work is 10 
bins dex~^. The number of particles in the inner most 
ellipsoidal shell at 0.1 Tg for this binning scheme ranges 
form ca. 2500 (7 = 0) to around 10^ (7 = 2). In the 
outer regions we have typically 0(10^) particles in the 
ellipsoidal shells. At a given resolution, the fluctuations 
are larger in regions with a flat profile (7 = 0) than in 
regions with a steep profile (7 = 2). They decrease as 
well in the outer regions where the profile is even steeper. 
Of course, it is expected to some degree that the shape 
finding algorithms will have difficulty in resolving the 
small density contrasts from shell to shell in a nearly 
homogeneous region (7 = 0), which explains the central 
fluctuations seen in this case. Using the r~^ weighting 
in methods S2 and E2 leads to a significant shift of the 
axis ratio towards higher values than expected. 

All methods find the correct orientation of the prin- 
cipal axes within the well resolved range. For example 
for method SI, the median deviation of |cos((5a) — 1|, 
where 5a is the angle between the measured and the cor- 
rect direction of the semi-major axis a, is O(10~^) for all 
three different profile types. For the other methods, the 
alignment is of comparable quality. 

4.3. Changing axis ratios - aligned orientation 

Of course, real halos do not have a constant axis ratio 
as a function of distance. Therefore, we varied the axis 
ratios according to the simple parametrization given in 
Equation (12). The axis ratios were fixed at rg to b/a = 
0.8 and c/a = 0.6 and the slopes of Sb/a = ^0.1 dcx^^ 



and 



c/a 



= —0.15 dex were used. The condition 



h/a > c/a was assured by capping the parametrization 
with minima and maxima. This is not a problem within 
our range of interest between 0.1-10 rg. The orientation 
of the principal axes is kept aligned with distance. 

Figure 2 shows again only the axis ratio b/a as a func- 
tion of distance. Methods SI and S3 still give the best 
results. The weighting by r~^ introduces a bias towards 
higher values. Now, the methods using an enclosed vol- 
ume (El and E3) start to show deviations as well. This 
is due to the enclosed integration volume picking up in- 
formation from inner regions of the halo, which has a 
different shape. This leads to a lag in distance until the 
axis ratios can adapt. For example, these deviations for 
methods El and E3 become larger if we choose the axis 
ratio to change faster as a function of distance, e.g. as 
for the axis ratio c/a with s^/a — —0.15 dex~ . 

In the case shown in Figure 2, the axis ratio decreases 
with distance which leads to too high values for methods 
El and E3. If we choose the axis ratio to increase with 




Figure 2. Measured axis ratio b/a as a function of distance for 
halos with different inner slope 7 = (top panel), 1 (middle panel) 
and 2 (bottom panel). The halos were initialized with a changing 
axis ratio b/a as a function of distance (thin dashed line). We fixed 
b/a = 0.8 at rs and used a slope of Si,/a = —0.1 dex~^. Methods 
SI and S3 give the best results, whereas methods El and E3 start 
to show systematic deviations. Methods S2 and E2 give again too 
high axis ratios. 



distance, then the methods El and E3 are giving too low 
values. 

Even the methods SI and S3 do not perfectly reproduce 
the expected values. Similar as in the case for methods 
El and E3, they lie above/below the expected value if the 
slope of the axis ratio is an decreasing/increasing func- 
tion of distance. The deviations for methods SI and S3 
are smaller than for methods El and E3. These system- 
atic deviations seen for methods SI and S3 are mainly 
due to the local mass density profile. In regions with a 
flat local proflle, the systematic offset is bigger than in 
regions with a steep mass density profile. Increasing the 
resolution and using a finer binning (i.e. smaller aver- 
aging volume) only marginally decreases the offset. For 
regions with a local mass density slope 7 ~ 1 — 2, the 
systematic deviations in the case of a varying axis ratio 
are of the order of 0(0.01) for axis ratios for methods SI 
and S3. 

Again, all methods find the correct orientation of the 
principal axes. The directional deviations are very small 
and similar to what we found in Section 4.2. 

4.4. Changing axis ratios - changing orientation 

In real halos, the orientation of the principal axis 
can change as a function of distance as well. This is 
parametrized again by using the functional form of Equa- 
tion (12). The axis ratios are kept changing as in Section 
4.3. Additionally, we vary the alignment of the princi- 
pal axes by setting {61,92,03) = (0.375,0.125,0.25) r at 
Vs, with T = 27r. For the slopes we use {se^jSe^^se^) — 
(0.05,0.05,0.05) T dex-\ 

Figure 3 shows the axis ratio b/a as a function of dis- 
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Figure 3. Measured axis ratio hja as a function of distance for 
halos with different inner slope 7 = (top panel), 1 (middle panel) 
and 2 (bottom panel). The halos were initialized with a changing 
axis ratio hja as in Section 4.3 (thin dashed line) but we varied 
the orientation of the principal axes as a function of distance as 
described in the main text. Methods SI and S3 give again the best 
results, whereas the other methods show systematic deviations. 
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Figure 4. Cosine of the the alignment angle 5a, the angle between 
the measured and the correct direction of the semi-major axis a, 
as a function of distance for halos of Figure 3. The methods that 
use an ellipsoidal shell as an integration volume are much better 
in recovering the local orientation of the matter distribution. 



tance. If the local mass density profile is well resolved, 
then methods SI and S3 are closest to the correct axis 
ratios. For the other methods we see some systematic 
deviations which depend on the details of the axis twist. 



The findings about the axis ratios are reflected as well 
in the orientation of the principal axes. In Figure 4, we 
show cos(5a) as a function of distance, where 8a is the 
angle between the measured and the correct direction of 
the semi-major axis a. All methods that use the enclosed 
integration volume show larger deviations in the orienta- 
tion than the methods using a homoeoid as integration 
volume. The best method for recovering the local orien- 
tation in the well resolved region in this case is 82 tightly 
followed by SI and S3. The deviations for the semi-major 
axis a, (S^,, are the largest. The deviations are smallest 
for the semi-minor axis c, i.e. we have shown the worst 
case in Figure 4. 

4.5. First conclusions 

We have experimented with many more mass density, 
shape and orientation profiles as well as different reso- 
lutions than shown here. The findings are always the 
same: using an ellipsoidal shell as an integration volume 
without or with r^j weighting (methods SI and S3) gives 
results that are closest to the expected value under con- 
trolled conditions in regions where the mass distribution 
is well resolved and the density contrast is high enough 
(i.e. no flat mass density profiles). 

Methods SI and S3 agree, since the weighting by r~^ in 
each shell is like dividing by a different constant in each 
shell, which does not affect the axis ratios. The absolute 
values of the eigenvalues of the shape tensor for method 
S3 change of course. Hence, our preferred method is the 
pure form without any weighting, i.e. method SI. All 
other methods lead to significant deviations that in de- 
tail depend on the mass density, shape and orientation 
profile. This makes it also impossible to come up with a 
correction scheme that works in all cases that would al- 
low to convert the measured axis ratios between different 
methods. 

5. HALOS FROM COSMOLOGICAL STRUCTURE 
FORMATION SIMULATIONS 

Now we turn to a study of halos in cosmological struc- 
ture formation simulations. In these halos, in addition 
to the change of the axis ratios and the orientation of 
the principal axes as a function of distance, we also have 
subhalos. 

The data are from a cosmological structure formation 
simulation, where we simulated several objects that will 
end up as Milky Way-sized objects at redshift z — Q. The 
simulations were run with the latest version of the gas dy- 
namics and TV-body adaptive refinement tree (ART) code 
(Kravtsov ct al. 1997; Kravtsov 1999; Kravtsov ct al. 
2002; Rudd ct al. 2008). ART includes 3-dimensional 
radiative transfer of ultraviolet (UV) radiation from in- 
dividual stellar particles using the optically thin variable 
Eddington tensor (OTVET) approximation (Gnedin & 
Abel 2001). It includes a non-equilibrium chemical net- 
work of hydrogen (Hi, Hii and H2) and helium (Hei, 
Hen and Hem) as well as non-equilibrium cooling and 
heating rates, which use the local abundances of atomic, 
molecular and ionic species as well as the local UV in- 
tensity (Gnedin & Kravtsov 2011). All these proper- 
ties are followed self-consistently during the course of a 
simulation. An empirical model for the formation and 
shielding of molecular hydrogen on the interstellar dust 
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Figure 5. Measured shape of the total matter distribution of a massive halo at z isi 2 in our cosmological simulation. In the top row 
wo plot b/a and in the bottom row c/a as a function of distance. In the left column we show methods S1-S3 (integration over ellipsoidal 
shell), in the right column methods E1-E3 (integration over enclosed ellipsoidal volume). Cases where we removed the subhalos are marked 
with a -, cases where they remained by a +. It is evident, that it is essential to remove the subhalos in order to calculate the local shape 
correctly. 



allows for more realistic star formation recipes based 
on the local density of molecular hydrogen (Gncdin & 
Kravtsov 2011). Also included in ART is metal enrich- 
ment and thermal feedback due to the Type II and Type 
la supernovae (Kravtsov 2003) as well as stellar feedback 
(Kravtsov & Gncdin 2005). Here, we use data at z ~ 2 
from a simulation that includes cooling and star forma- 
tion (simulation series A). Further details are presented 
in an accompanying paper (Zcmp ct al. 2011). 

Figure 5 shows the shape of the total matter distribu- 
tion of a massive halo at z « 2. The distance is nor- 
malized by r2oob = 101 kpc, the radius that encloses a 
spherical volume such that the average enclosed density 
is 200 times the background density at that epoch. The 
halo has a total mass M2oob — 1-13 x 10^^ M0 and con- 
tains 6.71 X 10^ gas volume elements, 5.38 x 10^ dark 
matter and 1.39 x 10^ star particles within r2oob- AH 
variants of the methods are shown with and without the 
subhalos from the resolution scale (0.003 r2oob) up to 2 

?'200b- 

Subhalos are removed by cutting out a spherical hole 
around the subhalo center with radius rti-unc- The spher- 



ical mass density profile of subhalos typically shows an 
uprise at large distances from their center due to the host 
halo. The location where the minimum mass density is 
reached defines the truncation radius rtj-unc- We investi- 
gated under controlled conditions the effects of cutting 
out holes of typical sizes of massive subhalos at different 
distances from the host halo center. The deviations for 
the measured axis ratios at the location of the subhalo 
can be a few percent for methods SI and El when com- 
pared to the smooth case. Alternatively, one could only 
remove particles bound to subhalos (Lau et al. 2011). 

For the halo shown in Figure 5, the most massive sub- 
halo has a mass of 9.02 x 10^ M0, rtrunc = 2.77 kpc and 
is located at a distance of 10.5 kpc w 0.1 r2oob from the 
host halo center. The total mass in all subhalos in this 
case is 3.71% and most of the subhalos are located in the 
outer region of the halo. 

Generally, the presence of massive subhalos leads to 
spikes in the axis ratios b/a and c/a when using an el- 
lipsoidal shell as an integration volume (methods S1-S3). 
The subhalos bias the measured axis ratios drastically 
at locations where they constitute a significant fraction 
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Figure 6. Measured median axis ratios 6/a (top panel) and cja 
(bottom panel) as a function of distance for the total matter dis- 
tribution of the 16 halos at 2 Ri 2 for all the methods without the 
subhalos. 



of the total mass in the elhpsoidal shell. These spikes 
are visible for all weight functions - most pronounced 
if no weighting or w{r) — r^^ is used and least pro- 
nounced for w{r) = r~^. This effect is still present, 
though weaker, when integrating over the enclosed ellip- 
soidal volume without any weighting (method El). Often 
it is claimed in the literature that using the weights r~^ 
or r~y^ in the shape tensor reduces the influence of sub- 
halos on the shape determination. This is true only if 
an enclosed integration volume is used (methods E2 and 
E3). 

Again, there is nearly no difference between methods 
S1-, S1-I-, S3- and S3-)- if there are only few or no subha- 
los present at that distance (i.e. in the inner region). If 
we integrate over the whole enclosed ellipsoidal volume, 
then the inclusion of the weighting by r~^ smoothes out 
the detailed shape features. Worse is using the weight 
r~^ which leads again to a systematic shift of axis ra- 
tios towards larger values in our case in addition to the 
smoothing already observed for the r^j^ weighting, i.e. 
shapes are determined as rounder than they actually are. 

The median shape of the 16 most massive halos at 
2; « 2 for all methods without the subhalos is shown in 
Figure 6. Taking the median is motivated by the similar- 
ity of our selected halos (within a factor of 10 in mass, see 
also Zemp et al. 2011). As already observed before, the 
methods SI- and S3- are nearly identical. The methods 
where we integrate over an enclosed ellipsoidal volume 
(E1-E3) are naturally smoother than when integrating 
over an ellipsoidal shell volume (S1-S3). As a conse- 
quence, the local shapes do not react as fast to shape 
changes in distance as in methods S1-S3, as seen for ex- 
ample for cj a. There is a lag in distance when compared 
with shapes determined by methods that use an ellip- 
soidal shell as integration volume. This is also visible for 



our single halo in Figure 5. For methods S2 and E2, the 
bias towards rounder shapes can be around 0.1-0.3 for 
both axis ratios. 

6. DISCUSSION 

A widespread method used in the literature is method 
E3 (e.g. Dubinski & Carlbcrg 1991; Allgood et al. 2006; 
Kuhlcn et al. 2007). By using an enclosed integration 
volume, this method picks up information from the inner 
regions that can have different shapes and orientation. 
If one is interested in the local shape, then we find that 
method SI is clearly a better choice than method E3. 

Method El (e.g. Katz 1991) is doing relatively well 
compared to its differential version SI. This is due to the 
fact that the contribution in the shape tensor (Equation 
(5)) is dominated by particles or volume elements with 
the largest distance from the center. This method also 
shows systematic shifts (see for example Figure 6) and 
smoothing when compared to method SI. Therefore, the 
differential version SI should be preferred over the El 
method that uses the enclosed ellipsoidal volume. 

Unfortunately, methods SI or S3 are not yet in 
widespread use in the literature. Kazantzidis et al. (2004) 
used method S3 and also found that using the enclosed 
volume is sensitive to the distribution of particles in the 
enclosed region. Unfortunately, they did not present the 
details of the tests in their work. Debattista et al. (2008) 
and Lau et al. (2011) are also advocating method SI. 
While Debattista et al. (2008) do not further motivate 
their choice, Lau et al. (2011) found from visual com- 
parison that using a differential method in 2 dimensions 
gives reliable ellipsoidal fits to X-ray isophotes. 

7. SUMMARY 

We have critically examined different methods for de- 
termining the local shape of matter distributions as a 
function of distance. Using the weights r~^ or r~[j in the 
shape tensor (Equation (9)) does not cure the problem 
arising due to the presence of subhalos. In contrary, it 
can lead to a systematic bias for the measured axis ratios 
even in smooth cases (Section 4). We think it is better 
to remove the cause of the problem (i.e. the subhalos) 
than to fight the symptoms with weight factors that make 
the physical meaning of the shape tensor unclear. Also 
when integrating over the whole enclosed ellipsoidal vol- 
ume, features get smoothed out and shape changes are 
lagging behind in distance. 

Therefore, our recommended method for measuring lo- 
cal shapes is removing the subhalos, using ellipsoidal 
shells as the integration volume and determining the 
shape through an iteration method as described in Sec- 
tion 3 that uses the shape tensor as defined in Equation 
(5), i.e. without any weight factors. 

In some cases one is interested to characterize the 
shape of an object with just one number, i.e. one is not 
interested in the internal structure and the local shape 
as a function of distance. Also, if the object is not well 
resolved (typically less than 0(10'*) particles/volume el- 
ements), calculating the local shape can be problematic. 
As a good practice, we recommend to have at least a few 
thousand particles in a bin when using ellipsoidal shells 
as integration volume. Therefore, if the internal struc- 
ture is not of interest or can not be properly resolved. 
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we advocate method El since this shows the least bias 
among the tested methods that use the enclosed ellip- 
soidal volume. 
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